How gravity and size affect the acceleration statistics of bubbles in turbulence 
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We report results from the first systematic Lagrangian experimental investigation in the previ- 
on : ously unexplored regime of very light (air bubbles in water) and large particles (D/n >> 1) in 

turbulence. Using a traversing camera setup and particle tracking, we study the Lagrangian accel- 
eration statistics of ~ 3 mm diameter (D) bubbles in a water tunnel with nearly homogeneous and 
isotropic turbulence generated by an active-grid. The Reynolds number (Rex) is varied from 145 to 
230, resulting in size ratios, D/n in the range of 7.3-12.5, where n is the Kolmogorov length scale. 
'-^ ' The experiments reveal that gravity increases the acceleration variance and reduces the intermit- 

tency of the PDF in the vertical direction. Once the gravity offset is subtracted, the variances of 
both the horizontal and vertical acceleration components are about 5 ± 2 times larger than the one 
t-|-H ' measured in the same flow for fluid tracers. Moreover, for these light particles, the experimental ac- 

^3 , celeration PDF shows a substantial reduction in intermittency at growing size ratios, in contrast to 

O ■ neutrally buoyant or heavy particles. All these results are closely matched by numerical simulations 

* ^5 ' of finite-size bubbles with the Faxen corrections. 

Oh 



I. INTRODUCTION 



^v^j ' Suspensions of particulate materials, drops or bubbles carried by vigorously turbulent flows occur frequently both 
^ , in the realm of natural phenomena (e.g. cloud formation) and in industrial applications (e.g. combustion in engines). 
In order to quantify the statistical properties of such suspensions a prototype problem is often considered: the one of a 
dilute suspension of spherical particles in incompressible, statistically homogeneous and isotropic turbulence [H, H| - 
^""^ i In this simple form the problem is defined by a set of three dimensionless parameters [Re\,T, 3], respectively the 
Reynolds number based on the Taylor scale of the carrying flow, the particle to fluid mass density ratio (r = p p /pf) 
and the particle to dissipative- length ratio (3 = D/n). Interesting theoretical questions concern how far the particle 
statistical properties (e.g. moments, probability density functions (PDFs), correlations) of position, velocity and 
acceleration depart from the Lagrangian properties of the fluid. The goal is to understand how such observables vary 
. . ■ as a function of the control parameters. 

Numerical studies have attempted to see how closely the approximate equations of Lagrangian dynamics - which 
were known for a long time - are able to capture the dynamics. Since full numericsfejj;. physalis 5], front tracking 
Q) are too expensive for hig h ReA, most simulations use a point particle model [JOf, a ^ so known as the Maxey- 
Riley-Gatignol model [13, lll| In the real world, there are many situations where the particle size is larger than the 
Kolmogorov length scale of turbulence (S 3> 1). Therefore, a considerable body of recent work (jl, [H-[i3) nas been 
dedicated to the characterization of these so-called finite-sized particles. Numerical simulations with Faxen corrections 
to the point-particle approach [14| correctly ca ptu re two important features from the experimental data for neutrally 
buoyant and heavy finite-sized particles @, 0, 13 the acceleration PDF of the finite-size particles in general 

show less intermittency than those of fluid tracers; and their acceleration variance decreases with increasing size ratio. 
However, for light particles (r < 1), the Faxen corrected numerics [14j remarkably indicate an opposite trend for 
the acceleration statistics at growing H: an initial increase of acceleration variance and intermittency, followed by a 
decrease. These predictions for light particles are awaiting experimental verification. However, these experiments are 
highly challenging in terms of the infrastructure needed, techniques and analysis, but are of great importance for a 
fundamental understanding of particles in turbulence. 

In this paper we present the first systematic experimental investigation in the regime (r <C 1,3 > 1) of the 
parameter-space, i.e. very light and large particles (see Fig. [1] for a summary of all the presently available mea- 
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FIG. 1: Parameter space of the density ratio versus size ratio for particles in turbulence, from data available in literature: A - 
Voth etal. H, o - Qureshi et al. [HI], + Volk et al. [l§], o Brown et al. [IB], * Gibert et al. [l]|, □ Volk et al. [13], ■ Martinez et 
al. [2(J, • Present work. Majority of previous studies have focused on Y > 1, here we explore the f< 1 case. 



surements). For such an investigation we use air bubbles, which are dispersed in a turbulent water flow. We track 
these bubbles using a traversing camera system which can perform 2D recordings of the vertical and one horizontal 
component of the bubble trajectories. The experiments are compared to numerical simulations based on the particle 
Lagrangian equations with Faxen correction [14( ■ We also experimentally study the effects of gravity, as it could be 
important at lower Re a, and very few numerical studies [lH have taken gravity into account. 

The paper is organized as follows: in section [TTl we describe in detail the experimental setup (sec. Ill Af . the bubble 
deformability and size distribution (sec. Ill Bp . Next, we present results on the velocity statistics (sec. and the 
acceleration statistics (sec. II VI) . The effects of gravity and size on the acceleration statistics are discussed in sections 
(sec. IIV A|) and (sec. HVBj) . followed by the conclusions in section [V] 



II. EXPERIMENTS AND ANALYSIS 



A. Experimental setup 

In our experiments, air bubbles in water, T w 10 -3 , are dispersed in a turbulent flow in the 8m high Twente Water 
Tunnel (TWT) facility (see Fig.EJa)). Nearly homogeneous and isotropic turbulence is generated by the flow of water 
through an active grid [2JJ. An optically transparent measurement section of dimensions 2x0.45x0.45 m 3 is located 
downstream of the active grid. We recently reported results on ~ 0.3 mm sized bubbles ([13, H3) using the same 
experimental facility (see also 0, [24[). Here, one important modification is made: the position of the active grid is 
switched from the position on top of the measurement section to its bottom. Furthermore, the direction of the water 
flow through the active grid is now upwards. The bubbles are generated by blowing air through capillary islands 
(diameter 500 /im) placed below the active grid. The bubbles rise through the measurement section along with the 
imposed water flow and escape through an open vent on top of the water tunnel. A surfactant (Triton X-100) is 
added to the tap water to reduce the bubble deformability [25( (more details follow in section III Bl) . Re>, is varied 
from 145-230 by changing the mean flow speed of water in the tunnel. The flow characterization is done using a 
hot-film probe placed in the center of the measurement volume [20T ] . The two-dimensional Lagrangian particle tracking 
experiments are carried out using a high-speed camera (Photron SA1.1) at an image acquisition rate of 5000 frames 
per second (fps) with a resolution of 768 x 768 pixels. The camera is focused on a 1-2 cm thick plane in the middle 
of the measurement section and the illumination is provided from the opposite side by a halogen light source placed 
behind a diffusive plate. The camera is mounted on a traverse system (Aerotech L-ATS62150 linear stage) which 
enables precise movement in the vertical direction at preset velocities. The detection of the bubble centers in the 
images is a non-trivial task because the bubbles often overlap or go out of focus (see Fig. L^b)). However, the circular 
Hough transform method [27| is successfully used to detect more than 90% of the bubbles (which are in focus) in the 
images. A 2D Particle Tracking code is then used to obtain the bubble trajectories over time. The data processing 
approach is the same as in Martinez et al. [2(|. The fitting window lengths (N) of the trajectory smoothing for the 
micro-bubble experiments (see Ref. [20l ] for more details) and present experiments have been chosen to be consistent. 
The fitting window length is N — 15 for the present experiments at 5000 fps acquisition rate, and N — 30 for the 
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FIG. 2: The Twente Water Tunnel facility: vertical water tunnel with nearly homogeneous and isotropic turbulence generated 
by an active-grid. Bubbles are dispersed from below through capillary islands and the flow is in the upward direction. The 
camera moves upward at preset speeds, along with the bubbles, allowing the measurement of long-duration trajectories. 



micro-bubble experiments at 10,000 fps acquisition rate. 



B. Bubble deformability and size distribution 

A surfactant (Triton X-100) is added to the tap water to reduce the bubble deformability. The small amount of 
surfactant (< 1 ppm) used in the present experiments is much below the critical micelle concentration, so the change 
in flow properties is negligible [26| . Although the addition of surfactant reduces the surface tension by a few percent, 
it leads to another competing effect, namely the suppression of bubble coalescence at the source of injection. This 
second effect is more dominant and as a result we see a reduction of the bubble size (25j . 

The Fig. [3] below shows the difference between bubbles in tap water and in a surfactant solution. We clearly see a 
big difference in the shapes and sizes of the bubbles. The deformation is greatly reduced in the presence of surfactant, 
thanks to the smaller size. Nonetheless, in Fig. [3l the bubbles still appear to be deformed and anisotropic, but 
the effect is much less than without surfactants: The addition of surfactants leads to a reduced deformability of the 
bubbles, not to a perfect spherical shape; the principal- axis deformation ratios are less than 1/2. The images in Fig. [3] 
are instantaneous snapshots from the experiments. It must be noted that the deformation of real bubbles is strongly 
time dependent, and the addition of surfactant significantly reduces this deformation over time. Thus, the addition 
of the surfactant renders the images amenable to processing, which is otherwise a nearly impossible task. 

The bubble diameter has a weak dependence on Re>, as shown in Fig. 2) We select the peak values of the distribution 
as the characteristic bubble diameter; D — 3.15,2.90,2.70,2.55,2.50 mm (with absolute deviations ~ ±0.3 mm, see 
Fig.@| at Re\ = 145, 170, 195, 215, and 230, respectively. This variation in ReA changes the Kolmogorov length scale 
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FIG. 3: Left panel: Rising bubbles of diameters ~ 5mm in tap water. Right panel: Bubbles of diameters ~ 3mm in surfactant 
solution (with < 1 ppm Triton X-100). 
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FIG. 4: Distribution function of bubble diameters, D, at different ReA. 



(77), and the corresponding size ratios are S = 7.3, 8.9, 10.2, 11.2, and 12.5. 



III. RESULTS - VELOCITY STATISTICS 



In the present experiments, the bubble motion and liquid mean flow are both upwards. The terminal velocity of 
the ~ 3mm bubbles in the absence of a liquid mean flow (still liquid) is determined from our experiments to be 25 
cms -1 . In the presence of a liquid mean flow, the bubbles attain their terminal velocity and then feel the turbulence 
of the surrounding liquid flow. In Fig. [5] we show the actual mean rise velocities of the bubbles (measured from the 
present experiments) compared to the liquid mean flow speeds. The bubbles experience a so-called 'slip' and are able 
to rise with velocities higher than the liquid mean flow. We quantify this "mean bubble slip velocity" by defining it 
as the difference between the bubble rise velocity and the liquid mean flow velocity. The variation of this mean slip 
velocity with Re a is shown in Fig. [SJ We see a decrease of the mean bubble slip velocity with increase in Re^. This 
decreasing trend might be due to the decreasing bubble size at higher Re>- 

We now present results on the probability distribution functions (PDFs) of the finite-sized bubble velocity obtained 
from the present experiments. Fig. (SJa, b) shows the PDF of the horizontal (x) and vertical (y) components of 
the normalized bubble velocity at the different Rc>, (and S) covered in the present study. We observe that the 
x-component velocity distributions (Fig. EJa)) follow sub-Gaussian statistics (flatness values: 2.27 — 2.78) and the 
y-component velocity distributions (Fig.^b)) show a slightly higher intermittency compared to the Gaussian profile 
(flatness values: 2.89 — 3.77). In previous experiments on finite-sized neutrally buoyant particles in turbulent von 
Karman flows, Volk et al. [171 ] have obtained sub-Gaussian velocity distributions with flatness around 2.4—2.6, but in 
the past, Gaussian-type flatness values (2.8—3.2) have also been reported [3j. The reason for deviations of the velocity 
distributions from Gaussianity is still an open question. 
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FIG. 5: Finite-sized bubble slip velocity at different ReA 
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FIG. 6: Velocity PDFs of the finite-sized bubbles at different ReA (and E), (a) Horizontal component, v x (b) Vertical component, 
v y (c) flatness F(v) versus ReA. The solid line in (a) and (b) represents a Gaussian distribution fit. The velocity PDF in (a) 
is mostly sub-Gaussian and in (b) we see more intermittency than a Gaussian profile. 



IV. RESULTS - ACCELERATION STATISTICS 



A. Acceleration statistics: Gravity effect 



We first address the effect of gravity on the acceleration statistics. Since the buoyancy is proportional to bubble 
volume, while the laminar viscous drag grows with the linear size of the bubble, it is clear that for growing bubble 
sizes and fixed turbulence intensity (growing S), buoyancy at some point shall dominate. The opposite is true for 
the case of fixed bubble size but increasing Reynolds numbers (again growing S) in our experiments, where buoyancy 
loses its dominance at higher Re\. In summary, it is acceptable to neglect the buoyancy force only for small bubbles 
S < 1 [20], or for large ReA, as we find later. 

One may expect that the buoyancy force will produce asymmetry: on the vertical component statistics because 
buoyancy will add up to the upward acceleration events and will subtract from the downward events. However, we 
find that the asymmetry is almost negligible. First, we find that the mean value of vertical component of acceleration 
(y) is essentially zero, as it is for the horizontal one (x): (a x ) ~ (a y ) = ± 0.2<?, where g is the acceleration due 
to gravity. Second, we observe that the probability density function (PDF) shape is only very weakly asymmetric; 
we indeed estimate the skewness, S{a,i) = (of )/(af ) 3//2 , and find ^(a^)! < 0.01 and \S(a y )\ < 0.1. The skewness is 
comparable to the values found for S < 2 bubbles (from now on called microbubbles) studied in the same setup [20| . 
We can conclude that if any asymmetry is present it must be very weak. 

The buoyancy force, however, produces a robust anisotropy: different statistics for the vertical and the horizontal 
acceleration components. This influence of gravity is clearly visible in the second statistical moments of acceleration, 
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shown in Fig Ufa), where the variance (a 2 ) is plotted. We find that (a 2 ) » {a 2 x ) for all S > 7.3 bubbles, while from 
the same figure it is evident that microbubbles have much closer variance values for the three cartesian components. 
If one calls a' y the vertical acceleration component in absence of gravity and assumes that {a 2 ) = ((a' y — g) 2 }, one gets 
{a 12 ) = (a 2 ) — g 2 , since (a' y ) = because of isotropy. As can be seen from figure EJa), this produces a reasonable 
collapse of the x-y data (a 2 ) — g 2 ~ (a 2 )- We emphasize that the collapse observed for the second order moment is 
non-trivial. It means that the statistical effect of gravity seems to be additive on the vertical direction with no effect on 
the horizontal component (see supplementary information for more discussions on this). An immediate consequence 
is that the effect of hydrodynamic forces coupling different cartesian directions, as for example the lift force, turns 
out to be unimportant here. 
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FIG. 7: (a) The single-component acceleration variance normalized by the gravitational acceleration for large bubbles and for 
micro- bubbles [2(ij |. (b) The normalized acceleration variance versus the size ratio. Legend : The present large bubble results 
(•) are represented with errorbars, along with a gravity offset for the vertical component a y . Open squares □: point-particle 
DNS of bubbles [3] triangles A & A: bubbles from DNS with Faxen corrections [HI], plus + : single experimental point for 
bubbles diamonds & ♦: DNS of neutrally buoyant particles with Faxen corrections [IJ. 



As mentioned above, we have found that the acceleration variance (a 2 ) in the vertical direction is augmented by an 
offset (or a correction factor) that depends on g 2 . Although we have found that the offset of ~ g 2 seems to work well for 
the present finite-sized bubbles, we do not yet have a physical understanding of the origin of this correction. When one 
examines the equations governing the motion of bubbles in a turbulent flow (Maxey-Riley-Gatignol equations) [ToLITljl 
including gravity; the corresponding acceleration variance offset would then be ~ Ag 2 , which does not work for the 
present data. Clearly, there seems to be a complex interplay between gravity and inertia, and more experimental and 
numerical work is needed before we can arrive at solid conclusion on the gravity correction factor. These corrections 
will be explored in detail in future work. 

The finite-sized bubble acceleration PDFs and the flatness values at two selected cases of Re> are shown in Fig. [5] 
The quality of statistics in the present experimental data is seen in the insets of Fig. H{a,b), which demonstrate that 
the fourth-order moments are well converged. Hence, we can calculate the flatness directly from the distribution 
without resorting to a fitting procedure (e.g. as in Volk et al. [I3)- In Fig. HJc), we observe that the flatness values of 
a y at ReA < 195 are less than those of a x . At higher ReA (> 215), the flatness values of a x and a y become comparable, 
and the corresponding PDFs show a nice collapse in Fig. [5fb) . The decreased intermittency in the vertical component 
a y for the low ReA (seen in Fig.[5]Ja)) is due to gravity. Apparently, there are two regimes: at lower ReA, gravity has 
an effect on the acceleration statistics in the vertical direction, which is no longer the case at higher ReA- 



B. Acceleration Statistics: Size effect 



Once the statistical influence of buoyancy has been disentangled, we now study how the gravity-less acceleration 
variance changes at increasing H, in other words we study the purely hydrodynamic size-effect on particles which are 
lighter than the surrounding fluid. It is convenient to look at the relative change of the bubble acceleration variance 
with respect to fluid tracers {a 2 f)- For dij we use here the acceleration of the microbubbles as found in Ref. [2fj| . 
because given their small size, 5 < 2, they behave almost like Lagrangian tracers. Furthermore, the ReA numbers 
studied in the present experiments are very close to the ones analyzed in Ref. [2(| and the flow conditions are the 
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FIG. 8: Normalized bubble acceleration PDFs of a x and a y components at ReA = (a) 145 and (b) 230, which are respectively 
the smallest and the largest ReA considered in our experiments. As expected, at the lower ReA the tails of the normalized PDF 
(in other words the flatness) of the vertical component are reduced compared to the horizontal one, while at the higher ReA 
the anisotropy is negligible. The insets in (a) and (b) show the fourth-order convergence test, (c) The acceleration flatness for 
the present bubbles (with errorbars) versus ReA- 

same. Fig. [TJd reports such a normalized acceleration variance (af) / (af j) versus the size ratio S. We see that both 
acceleration components reach a level of 5 ± 2 times the variance measured in the same flow for fluid tracers. They 
are also in agreement with the previous single experimental datapoint of Ref. (l9| . To have an interpretation of 
this measurement we compare it with the results from numerical simulations at similar Re,\. Fig. [TJd reports results 
from two different types of Lagangian particle simulations, first the so called point-particle (PP) simulation that only 
takes into account the hydrodynamic effects of added mass and Stokes-drag, and a second simulation which adds on 
Faxen corrections (FC) to the mentioned terms [2 [HI- The PP model predicts for S — )• oo an asymptotic limit of 
normalized acceleration variance which is 9 times that for tracers (as a result of the dominance of the added-mass 
term). From Fig[7jD, we deduce that for 3 « 10 in the PP model, it is about 15 times the value for tracers, reflecting 
that we are not yet in the asymptotic limit. For larger 2 the results from the PP model indeed seem to approach the 
asymptotic value of normalized acceleration variance of 9. Faxen corrections to the added mass term reduce the value 
at H = 10 which is about 7. Fig[7jD clearly shows that numerical simulations of Faxen corrected finite-sized bubbles 
show a good agreement with the experimental measurements (contrary to the PP model which overestimates the 
result). The present results are the first systematic measurements to confirm that the acceleration variance increases 
with finite-size for bubbles, matching the FC numerical simulations [lij]. Note again that the trend is very different for 
the case of neutrally buoyant particles and heavy particles, for which the normalized acceleration variance is always 
< 1 0[H (Fig.Hb). 

Making statements on the fourth order statistical moments of acceleration based on experimental measurements is 
a delicate endeavor, see e.g. the detailed analysis by Volk et al. [lj|. To study the finite-size effects on the bubble 
acceleration, in the following we will focus only on the PDF shape of the normalized x— acceleration, a x j (a x ) Tlas , 
which is not directly affected by gravity. In Fig. HJa) first we plot such a curve for microbubbles [2(| and see that its 
shape falls on the one for fluid tracers, and on S < 2 bubbles, from DNS at similar ReA- This is an evidence of both 
the above mentioned passive nature of microbubbles and of the similarity between turbulence realized in the Twente 
Channel flow and the one produced in homogeneous and isotropic DNS. In the same panel we see that, in sharp 
contrast, the finite-sized bubbles show a strongly reduced intermittency. It is the first time that such a substantial 
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FIG. 9: The normalized acceleration PDFs of the horizontal component \a x \, from the present bubble experiments (EXP), 
microbubble measurements [2(j, and numerical simulations (DNS) with the Faxen model at Re^ = 180 (TJ. (a) Comparison 
with results for point-like particles; DNS Tracers, microbubbles, and Faxen DNS of bubbles for S < 4. (b) Comparison with 
Faxen DNS of larger bubbles, 3 > 8. 



change in intermittency at growing size ratios is experimentally observed: neither for solid neutrally buoyant particles 
[HIE]]) nor f° r heavier bubbles [l3[ was it detected before. 

In order to clarify further the magnitude of such an effect, in Fig.[9^b) we compare the acceleration PDF from the 
DNS simulations with Faxen corrections and the present experiments. There is an initial increase of flatness in the 
DNS simulations for S < 8, which probably reflects the limitations of the Faxen corrections. For 5 > 7 — 8 both 
the numerics and the experiments show a significant reduction of the tails of the PDF. However, the DNS appears 
to underestimate its functional behavior by approximately a factor 2-3 in the size parameter H. The reason for this 
discrepancy is presently unclear. One reason could be that the simulations neglect the two way coupling and are 
just approximated in the implementation of Faxen terms. Another possible reason for the discrepancy which deserves 
further study is the deformability of real bubbles. The deformation process absorbs/releases energy from/to the 
turbulent environment, a process which may have an effect on acceleration statistics. These issues motivate further 
investigations to better understand light particles in turbulence. 



V. CONCLUSION 



We performed measurements of Lagrangian acceleration in the previously unexplored regime of large (compared to 
rj) and very light (with respect to pf) particles in turbulence. Bubbles of size ^3 mm diameter (D) were tracked in 
turbulent flow conditions in a water tunnel. The explored range of Reynolds number (Re\) and size ratios, 3 = D/rj 
were 145-230, and 7.3-12.5, respectively. 

Gravity produces anisotropy in the acceleration statistics of the vertical component - it adds a g 2 offset to the 
variance, and decreases the intermittency of PDF, at lower Re>. It was found that the interaction between gravity 
and inertia is complex, and this deserves further study. 

The acceleration variances and the intermittency clearly indicate the finite-size effect, and the results are in good 
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agreement with DNS simulations with the Faxen corrections. To improve the current understanding, in the future we 
plan to study rigid hollow spheres in the turbulent water flow, and to vary their diameters, D, at fixed Re^ numbers. 
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